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The brightest blazars detected by the Fermi Gamma-ray Space Telescope Large Area Telescope 
(FermifLALT) show significant breaks in their spectra at a few GeV. The sharpness and the position 
of the breaks can be well reproduced by absorption of 7-rays via photon-photon pair production 
on He II and H I Lyman recombination continua (LyC) produced in the broad-line region (BLR). 
Using 138 weeks of LAT observations of the brightest GeV blazar 3C 454.3 we find a power- 
law dependence of the peak energy on flux and discover anti-correlation between flux and the 
column density of the Hell LyC which is responsible for absorption of the > 2.5 GeV photons 
in this object. The strength and the variability of the absorption implies the location of the 7-ray 
emitting zone close to the boundary of the high-ionization part of the BLR and moving away 
from the black hole when the flux increases. A combination of the GeV breaks with the detection 
of a few powerful blazars in the TeV range puts strong constraints on the BLR size. Additional 
spectral breaks at ^100 and ^400 GeV due to absorption by the Balmer and Paschen lines could 
be detected by the Cherenkov Telescope Array. 
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1. Fermi spectra of blazars 



The Fermi/LAI has detected a few hundred blazars in the 100 MeV-100 GeV range [Q]. The 
high photon statistics from the brightest sources allows also to study their spectra with unprece- 
dented accuracy. These studies have shown that the spectra of high-luminosity sources, flat-spectra 
radio quasars (FSRQs) and low-energy synchrotron peaked BL Lac objects, are much better de- 
scribed by a broken power law than by a simple powerlaw or any smoothly curved models ^ ^ 
(see Fig. |l]a). The works dedicated to 3C 454.3 [§, 0] found the spectral hardness-flux correlation 
at a nearly constant break energy. 

The observed spectral breaks are too sharp to be associated with the cooling or the Klein- 
Nishina effects [||]. External Compton scattering on a truncated electron spectrum ^ also cannot 
possibly reproduce the break sharpness. Naively, one would think that the break energies, mostly 
lying below 10 GeV (as measured in the object frame, see Fig. |I]a), are too small to be consis- 
tent with the y-ray absorption due to photon-photon (77) pair production in the broad-line region 
(BLR), as the strongest BLR line, Lya, absorbs 7-rays starting only from 25.6 GeV [^. However, 
the BLR contains not only hydrogen, but also e.g. helium, and emits not only lines, but also strong 
recombination continua. The observed break energies [||, [T8|] favour the He II and H I Lyman re- 
combination continua (LyC) as the main opacity source. A relatively high opacity in the He II LyC 
observed in a number of bright blazars implies the location of the 7-ray emitting region within the 
highly ionized inner part of the BLR [HI 
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Figure 1: (a) Spectral energy distribution of a few blazars as observed by FermilhP^ during first 6 months 
of its operation. The best-fit broken power law and a power law with the absorption by the He II and H LyC 



(double-absorption) are shown by the dashed and solid lines, respectively. From |18|. (b) Spectral energy 
distribution of 3C 454.3 at low and high fluxes (intervals B and H, see Fig. ^ as well as averaged over the 
whole observation period of 2.5 years. Here the underlying continuum is modeled as the lognormal function 
(shown by the dotted curve). The models modified by double-absorption are shown by solid curves, while 
the dashed curve is for absorption by BLR at a fixed ionization parameter log ^ = 2.5. From pQ]. 
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2. Broad-line region and GeV opacity 

It is well known that BLRs around quasars emit lines associated with different ionization 
stages, with the broader high-ionization lines being produced closer to the central black hole [p^]. 
Reverberation mapping also demonstrated a strong anti-correlation between the line width and the 



time delays of the line response to the continuum variations ||17|]. The BLR size, as measured 
from the C IV 1549A line delays, scales with the accretion luminosity as /? civ ~ O.IlIj"^ pc [p^. 
However, the delays in the high-ionization Hell 1640A line as observed in Seyferts are 3-5 times 



smaller [13, 16]. On the other hand, the B aimer lines give sizes 2-3 times larger. It is also worth 



mentioning that the scaling [ |12| | extrapolated to the high-luminosity quasars is based on one object 
(S5 0836-1-71) and one line (C IV 1549A), while the BLR size estimated in the Seyfert-luminosity 
range has more than an order of magnitude spread. Thus the actual BLR size in quasars is highly 
uncertain. The strongest BLR line, Lya, does not show any variability on the timescale of years in 



the four quasars in which it was observed [12], implying that it is produced much further away that 
C IV 1549A. This has important implications for the opacity of the GeV photons. 

Using XSTAR we generated [ p^ ] a grid of photoionization models of the BLR clouds assumed 
to be simple slabs of constant gas density nu, of fixed column density at A^h = 10^^ cm^^, and a 
clear view to the ionizing source. We varied the ionization parameter = L/{r^nii) from lO"'^ 
to 10^^. The resulting spectra are shown in Fig. ^a. BLR radiation consists mostly of lines and 
recombination continua, which in the UV range produce narrow, line-like features, because the 
temperature of the photoionized regions is typically much below the corresponding ionization po- 
tentials of H and He. 

If the 7-rays are produced within the BLR, they propagate through a roughly isotropic photon 
field. The absorption strength by the isotropic Une photons of energy Eq can be characterized by 



the Thomson optical depth QISO 



A 7 ^line n . , „ ^line,45 10 cV 
TT=NphOT^——^-—ROT^\lO—^ —, (2.1) 

where A/ph is the column density of line photons along the line of sight, Oj is the Thomson cross- 
section, Liine is the line luminosity, and R is the typical size. (Here we define Q = 10^2x in cgs 
units.) The energy-dependent opacity is 

Tyy{E ,Eo)= NpkOyyis) = Tt . (2.2) 

Here Oyy{s) is the angle-averaged 77 cross-section, which has a threshold ats = £'£'o/(meC^)^ = 1, 
i.e. 

£th = 261GeV/(£o[eV]), (2.3) 

rapidly grows to the maximum of '--'0.21 ax at the energy of '^3.5 times the threshold, and then 
slowly decreases roughly inversely with energy High-energy photons transmitted through 

the soft line photons are attenuated as oc exp ( — T^y (£ ,Eo)). Approximating the spectrum by a power 
law, we can estimate the break in the power-law index at Eth produced by an individual strong line 
(or recombination continuum) as 

dlnexp(-Tyy(£,£o)) Onis)/aj Tt 

Ai — — — Txmax — . {ZA) 

alnE Ins 4 
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Figure 2: (a) Spectra Fe — dN/d\nE of the BLR clouds of column density A^h = 10 cm^ and various 
ionization parameters ^ photoionized by a quasar continuum shown by a smooth curve and consisting of a 
multicolor disc with the maximum temperature 10^ K and a power law tail of photon index F = 2 extending 
to 100 keV with 10% of the total luminosity. We use photoionization code XSTAR vs. 2.2 The upper 



axis show the threshold energy for 77 reaction corresponding to the BLR photon, given by Eq. (2.3). (b) 
Angle-average cross section (in units of Thomson cross section) for y-rays due to 77 absorption on the BLR 
photons with the spectra shown in panel (a). Jumps in opacity are clearly seen at energies corresponding to 
the strongest lines and recombination continua. The 0.3-0.5 GeV jump for large ^ is produced by O VIII 
16-19AA lines, the jump at low £, is due to the O VII 22A line complex. The jump at ~5 GeV visible at all 
ionizations is due to the He II LyC and Ly lines at 40-60 e V. The 20-30 GeV jump seen at low ^ is produced 
by H I Ly lines and continuum and the C IV 1549 A line. From 



Equation (2A) shows that in powerful quasars with the Lya luminosity of about 10^^ erg s^^, 
the 7-rays above 30 GeV have troubles of escaping unless the density of the line photons is suffi- 
ciently small (e.g. if they are produces in a large volume of /? > 3 pc), or if the 7-rays are produced 
outside the BLR. As the BLR size R (in a given line) scales with the total luminosity L^/^ [|^, the 
opacity also scales similarly (assuming a constant covering factor of BLR clouds) 

(2.5) 

and obviously the 77 absorption can be important only in the high-luminosity sources such as 



FSRQs. Another important consequence of Eq. (2.1) is that it is not the luminosity, but L/R ratio 



(i.e. compactness) that determines the role of a specific line in absorption of the 7-rays. 

Because of strong stratification of the BLR extending over two orders of magnitude in distance 



[14, hg], any model assuming that all lines are produced at the same distance from the central 



source [ p^ , [T9| , 22] strongly underestimates the 77-opacity by the relatively weak, high-ionization 
lines produced in a small volume (if the 7-rays are produced closer to the centre) and overestimates 
the opacity for multi-GeV photons from the strong, low-ionization lines which are produced in a 
larger volume. 
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The total 77-opacity associated with a given BLR shell at a fixed ionization can be determined 
by a convolution of the opacity for the line radiation with the BLR photon column density spectral 
distribution 

Zyy{E) = N^^Oyy{E) = Tt = / Oyy{s) -T^dE^. (2.6) 

The 77 cross sections weighted with the BLR spectra from Fig. ^ are shown in Fig. If ion- 
ization is high, the dominant photon source is the He II complex at 40-60 eV, which produces a 
break in opacity at 4-7 GeV. Above 10 GeV, the opacity is a rather smooth function of energy with 
additional absorption coming from He I, hydrogen Ly lines and recombination continuum, and C V 
2274A. In the 0.3-0.7 GeV region the opacity has another break due to high-ionization lines of 

VII. In the low-ionization environment, the strongest absorption is produced by hydrogen Ly 
lines and LyC, forcing a break in opacity at about 20-30 GeV. At the same time, the break due to 
helium is also clearly visible at a few GeV. The opacity is nearly flat up to 1 TeV because of the 
contribution from additional lines: Mgll 2800A, Ha, and He I 10832A. As we argued above, the 
complete model for the opacity has to be necessarily multi-zone and include radial stratification of 
the BLR and probably the angular dependence of the BLR radiation, which can affect the threshold 
energies. 

3. Blazars and absorption by He ii and H i LyC 

For the first analysis [|l8|], we have selected several brightest FSRQs from the sample of 12 
objects in Table 1 of and have chosen the same 180-day interval for easier comparison. We 
have developed our own software for data analysis and compared its results for the simple broken 
power law models with those obtained by the standard maximum likelihood analysis tool GTLIKE 
in For a few objects the fits with the power-law are acceptable, while other show clear breaks 
(see Fig. [T]a) that are rather well described by the broken power law model (see details in [|T8|]). 

To check the hypothesis that the BLR photons are responsible for the breaks, we have applied a 
simplified model of the BLR consisting of two strongest features: He II LyC at 54.4 eV and H I LyC 
at 13.6 eV. We call this double-absorber (DA) model. The free parameters are the optical depths 
iHe and Th in these"lines". The observed absorption threshold were redshifted by the appropriate 

1 + z factor. Assuming that the underlying continuum is a power-law, we obtained a good fit (see 
Fig. 0a). In RGB J0920-I-446 (with z. = 2.19), the absorption is seen only at high energies with 
the break energy corresponding to the pair-production threshold on H I LyC, 19.2 GeV/(l +z) ~ 6 
GeV, while in other cases the first break is close to the threshold for the absorption on He II LyC, 
4.8 GeV/(l and there are indication of the second break at higher energy due to Hi LyC. 

The FSRQ 3C 454.3 at redshift z = 0.859 is by far the brightest blazar in the 7-ray range 
during the lifetime of Fermi. The source shows variations in brightness by two orders of magnitude 
(see Fig. ^. LAT has detected more than hundred thousand photons from this source. With such 
a rich photon statistics, it is worth trying to study the spectral variation at various flux levels, to 
constrain the shape of the underlying spectrum, and to search for correlations between flux and the 
77 absorption depth. The last aim is the most important one, as it allows to learn about variations 
in the location of the 7-ray emitting region. 
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Figure 3: The light curve of FSRQ 3C 454.3 above 1 GeV as observed by Fermi/LAT. Adapted from 

We have split the data into the ten time intervals covering periods of more or less constant 
flux levels and analyzed the spectrum in each interval. The spectra for the high- (interval H) and 
low-flux (B) intervals as well as the average spectrum are shown in Fig. [l]b. Clear deviations from 
a simple power-law spectrum are obvious, and addition of the DA still does not provide a good fit 
at all flux levels. Thus instead of the power-law to describe the underlying continuum we have tried 
a lognormal function EFe °<: 10~'°^ (^/ViA)/°'in plus the DA. The fits with this model are superior 
compared to the previously used simple phenomenological models, especially at high fluxes with 
high photon statistics. The data can be also well fitted with 77 absorption produced by the full 
BLR spectrum of high-ionization (see dashed line in Fig. |l]b). A large contribution to the of the 
spectral fits comes actually from the high-energy tail. This might be an evidence for variations of 
the absorber optical depth during the observation and the fact that the sum of the spectra with low 
and high opacities differs from the spectrum at a fixed intermediate opacity. A probable excess of 
photons above 15 GeV (especially for the average spectrum) can also indicate that the spectrum is 
a superposition of emission states with different opacities due to the extended or moving emission 
zone. This proposal is supported by the spectral variability seen during the latest flare in 2010 
November [^, where photons above 10 GeV arrive mostly at the end of the flare, clearly indicating 
the position of the 7-ray emitting region further away from the black hole. 

3C 454.3 demonstrates a highly significant correlation of the peak energy Speak with flux (sim- 
ilar to the hardness-flux correlation [^), which is well represented by a power-law (Fig. ^): 

Speak = (46 ± 2) MeV f Oi,o±o.o4^ (3_1) 

where F = EFe at 1 GeV. The opacity Th ~ 3-7 did not vary significantly, but we observed an 
anti-correlation between Tne and the flux (see Fig. ^): 

THe = (2.61 ±0.22) -(1.15 ±0.47) logF^g. (3.2) 
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Figure 4: (a) Dependence of the peak (in EFe) of the best-fitting lognormal (with a\n = 2.7) + DA model 
on flux at 1 GeV as observed by Fermi/LALT from 3C 454.3. (b) Dependence of the Thomson optical depth 
in the He II LyC (defined by Equation (12. Ih) on flux at 1 GeV. The dashed line shows the best-fitting relation 



given by Equation (3.2). Adapted from [M 



The strength of absorption at GeV energies together with the estimations of the luminosities 



in LyC in principle constraints the BLR size (from Eq. p.l| ), where most of this soft radiation is 
produced. Thus, we get: 

/?He~0.22LHe,44(THe/3)pC, /?H ~ 3.7 Lh,45 ( W7) PC (3.3) 

The luminosities of 3C 454.3 in Lya and Hell 1640A are lO'*^ erg s"^ and 6 x 10"^^ erg s"^ 
respectively [^]. Assuming equal luminosities in LyC and Lya, the resulting Ru is much larger 
than the usually quoted sub-pc size of BLR. It is, however, consistent with the absence of Lya 



variability in powerful quasars [12|. The Hell LyC luminosity is probably larger than 10"*^ erg s^\ 
but is model dependent and therefore the estimate of /?He is less certain. 

A rather large ratio Thb/th (between about 1/4 and 1) and the observed Tug-F anti-correlation 
indicate that the /-ray emission region lies close to the boundary 7?He of the complete ionization of 
helium and moves out at higher luminosity. A small distance to the beginning of the /-ray emitting 
region alleviates the problems created by observations of the fast variability from 3C 454.3 [23, 8|]. 



4. Absorption features in the CTA range 

As we have shown in Fig. |2[ the strongest spectral features responsible for the 77 absorption 
and the breaks in the Fermi/LAI band are LyC of hydrogen and ionized helium. At energies above 
30 GeV, there are no significant spectral features expected, just because there are many closely- 
spaced BLR lines that produce rather smooth absorption profile. However, these simulations only 
covered a rather narrow set of parameters. The BLR not only has a large spread in the ionization 
parameter, but the cloud density may also vary by orders of magnitude. For simplicity we fixed 
logi§ and looked at the effects of the varying density (see Fig. ^. We see that in the high-density 
regime, the BLR spectrum below a few eV is dominated by lines of the hydrogen series. Strong 
Balmer lines (Ha, Hj8) give a rise to the 77 opacity at 100-140 GeV and the Paa line is responsible 
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Figure 5: (a) Same as Fig. |^a, but here we fix the ionization parameter at log^ = 2.5 and vary the density 
of the clouds from «h = 10"^ (at the top) to 10^^* cm^^* (at the bottom) with steps of a factor of 10. (b) 
Corresponding 77 cross sections weighted with the BLR spectrum are shown by the solid, dotted, dashed 
and dot-dashed curves, respectively. In the high-density regime there are additional jumps in opacity at ~100 
GeV due to the Balmer Unes and at ~400 GeV due to the Paa line. 

for another rise at ~400 GeV. Thus spectral breaks at these energies are expected if multi-GeV pho- 
tons are transmitted through such a BLR. Detection of such breaks with the Cherenkov Telescope 
Array (CTA) will give a clue to the structure and physical conditions within the BLR. 



5. Conclusions 

We have shown that the spectral breaks detected by the Fermi/LAI in powerful blazars at a 
few GeV can be naturally explained by the 77 absorption on the hydrogen and He II LyC. The 
strength of the absorption ultimately proves that the blazar zone has to be located within or close to 
the high-ionization zone of the BLR. Variability in the column density of He II LyC that produces 
absorption above ~3 GeV in the brightest GeV blazar 3C 454.3 implies that the 7-ray emitting 
region is located close to the boundary of fully ionized He zone and is moving away from the black 
hole when the flux increases. The strength of the Lya in this object allows to estimate the BLR size 
in LyC, which seems to exceed significantly the estimations previously discussed in the literature. 
A relatively large BLR size dilutes the density of the BLR photons and allows the multi-GeV 
photons to escape. Detections of a few powerful blazars in the TeV range thus does not necessarily 
means that the 7-ray emitting region is located outside the BLR as argued recently [Q, 21, ^]. On 
the contrary, the GeV breaks are impossible to produce far away from the BLR as there are no 
enough photons available for absorption and there is no any other physical mechanism that can 
produce sharp spectral breaks. High-density clouds embedded within the BLR can produce strong 
lines of hydrogen series which can absorb multi-GeV photons. Detection of the spectral breaks at 
~100 and ~400 GeV with the CTA will put interesting constraints on the BLR physical conditions 
and the location of the 7-ray emission zone. 



8 



Fermi Observations ofBlazars 



Juri Poutanen 



Acknowledgments 

This research was supported by the Academy of Finland grant 127512 and the Magnus Ehrn- 
rooth foundation. The research made use of public data obtained from the Fermi Science Support 
Center. We thank Maxim Barkov, Markus Boettcher, Evgeny Derishev, Chuck Dermer, Gabriele 
Ghisellini, Julian Krolik, Amir Levinson, Fabricio Tavecchio, and Dmitry Yakovlev for useful dis- 
cussions and Tim Kallman for his help with XSTAR. 

References 

[1] A. A. Abdo et al., Early Fermi Gamma-ray Space Telescope observations of the quasar 3C 454.3, ApJ 
699 (2009) 817 

[2] A. A. Abdo et al., Bright active galactic nuclei source list from the first three months of the Fermi 
Large Area Telescope all-sky survey, ApJ 700 (2009) 597 

[3] A. A. Abdo et al., Spectral properties of bright Fermi-detected blazars in the gamma-ray band, ApJ 
710 (2010) 1271 

[4] A. A. Abdo et al., Fermi Gamma-ray Space Telescope observations of the gamma-ray outburst from 
3C454.3 in November 2010, ApJ 733 (201 1) L26 

[5] M. Ackermann et al., Fermi Gamma-ray Space Telescope observations of gamma-ray outbursts from 
3C 454.3 in 2009 December and 2010 April, ApJ 721 (2010) 1383 

[6] F. A. Aharonian, Very high energy cosmic gamma radiation: a crucial window on the extreme 
Universe, World Scientific Publishing, River Edge, NJ, 2004 

[7] J. Aleksic et al., MAGIC Discovery of Very High Energy Emission from the FSRQ PKS 1222+21, ApJ 
730 (2011) L8 

[8] G. BonnoU et al.. The j-ray brightest days of the blazar 3C 454.3, MNRAS 410 (201 1) 368 

[9] R. J. Gould & G. P. Schreder, Pair production in photon-photon collisions, Phys. Rev. 155 (1967) 1404 

[10] H. T. Liu & J. M. Bai, Absorption of 10-200 GeV gamma rays by radiation from broad-line regions in 
blazars, ApJ 653 (2006) 1089 

[11] T. Kallman & M. Bautista, Photoionization and high-density gas, ApJS 133 (2001) 221 

[12] S. Kaspi et al.. Reverberation Mapping of High-Luminosity Quasars: First Results, ApJ 659 (2007) 
997 

[13] K. T. Korista et al.. Steps toward determination of the size and structure of the broad-line region in 
active galatic nuclei. 8: an intensive HST, lUE, and ground-based study ofNGC 5548, ApJS 97 
(1995) 285 

[14] J. H. Krolik, Active galactic nuclei: from the central black hole to the galactic environment, Princeton 
University Press, Princeton, NJ, 1999 

[15J D. E. Osterbrock, G. J. Ferland, Astrophysics of gaseous nebulae and active galactic nuclei. 
University Science Books, Sausalito, CA, 2006 

[16] B. M. Peterson & A. Wandel, Keplerian Motion of Broad-Line Region Gas as Evidence for 
Supermassive Black Holes in Active Galactic Nuclei, ApJ 521 (1999) L95 



9 



Fermi Observations ofBlazars 



Juri Poutanen 



[17] B. M. Peterson & A. Wandel, Evidence for Supermassive Black Holes in Active Galactic Nuclei from 
Emission-Line Reverberation, ApJ 540 (2000) L13 

[18] J. Poutanen & B. E. Stem, GeV breaks in blazars as a result of gamma-ray absorption within the 
broad-line region, ApJ 717 (2010) L118 

[19] A. Reimer, The redshift dependence of gamma-ray absorption in the environments of strong-line 
AGNs, ApJ 665 (2007) 1023 

[20] B. E. Stern & J. Poutanen, Variation of the JJ opacity by the He II Lyman continuum constrains the 
location of the y-ray emission region in the blazar 3C 454.3, MNRAS, in press (201 1) 
[arXiv:1105.2762] 

[21] Y. T. Tanaka et al., Fermi Large Area Telescope Detection of Bright J-Ray Outbursts from the Peculiar 
Quasar 4C -\-2L35, ApJ 733 (201 1) 19 

[22] F. Tavecchio & D. Mazin, Intrinsic absorption in 3C 279 at GeV-TeV energies and consequences for 
estimates of the extragalactic background light, MNRAS 392 (2009) L40 

[23] F. Tavecchio, G. Ghisellini, G. Bonnoli, G. Ghirlanda, Constraining the location of the emitting region 
in Fermi blazars through rapid y-ray variability, MNRAS 405 (2010) L94 

[24] F. Tavecchio et al.. On the origin of the gamma-ray emission from the flaring blazar PKS 1222+216, 
A&A, in press (201 1) [arXiv :1104.004 8] 

[25] B. J. Wills Hubble Space Telescope Sample of Radio-loud Quasars: Ultraviolet Spectra of the First 
31 Quasars, ApJ AAl (1995) 139 

[26] A. A. Zdziarski, Saturated pair-photon cascades on isotropic background photons, ApJ 335 (1988) 
786 



10 



